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Abstract 

We consider the fermion mass spectrum in the strong coupling vortex phase (VXS) of a lattice 
fermion-scalar model with a global U{1)l x U{1)r, in two dimensions, in the context of a recently 
proposed two-cutoff lattice formulation. The fermion doublers are easily made massive by a strong 
Wilson-Yukawa coupling, but in contrast with the standard formulation of these models, in which 
the light fermion spectrum was found to be vector-like, we find massless fermions with chiral 
quantum numbers at finite lattice spacing. When the global symmetry is gauged, this model is 
expected to give rise to a lattice chiral gauge theory. 


1 Introduction 


Wilson-Yukawa models were extensively studied some years ago as a possible lattice formulation of 
chiral gauge theories |^. The idea was to decouple the lattice doublers from the light spectrum by 
introducing the so-called Wilson-Yukawa term, which contains a charged scalar held in a way that it 
is manifestly gauge invariant under a chiral gauge transformation. With the addition of this “Higgs” 
particle, it was expected that doublers would decouple, while chiral fermions would remain in the light 
spectrum. This expectation was based on the numerical hnding of a strong Yukawa-coupling regime, 
both in the symmetric (PMS) and broken (FMS) phases, where the fermion masses are not proportional 
to the vacuum expectation value of the scalar held, as in the perturbative regime, and doublers can get 
masses of the order of the cutoh. However, for the same reason, also the light fermions are massive in 
this phase (although their mass could be tuned to zero by tuning the bare mass). More problematic was 
the fact that the light fermions were found to have only vector-like couplings to the gauge helds. The 
physical picture that emerged from these results was that in the strong coupling phase, composites of 
the light fermions and the scalar particles were formed, generating effective mirror states and rendering 
the theory vector-like. This situation was found already in the quenched approximation. 

It is important to realize that the failure of Wilson-Yukawa models is also very relevant for all 
those lattice proposals of chiral gauge theories that explicitly break the gauge symmetry. The reason 
is that, once gauge invariance has been broken for a compact group, the gauge non-invariant theory 
is exactly equivalent to a gauge-invariant theory with new charged scalar degrees of freedom |^. The 
scalar degrees of freedom are nothing more than the pure gauge transformations. The effective action 
of these scalars differs in general from one proposal to another, but they are coupled as long as gauge 
invariance is broken in the original formulation and so they must be dealt with. In particular, the 
simplest formulation of chiral gauge theories, which introduces a standard Wilson term to decouple 
doublers, breaks gauge invariance and is in one-to-one correspondence to the Wilson-Yukawa models 
considered in |l|. Equivalently all proposals that break gauge invariance must have a corresponding 
Wilson-Yukawa picture. It is important to understand to what extent the failures of these models in 
providing a regularization for chiral gauge theories are circumvented by other proposals. 

Recently a new method has been proposed to regulate chiral gauge theories on the lattice [^. In 
this approach, the fermions live on a lattice of spacing / and are coupled to gauge link variables that 
are constructed by an appropriate smooth and gauge-invariant interpolation of gauge conhgurations 
0 that are generated on a coarser lattice of spacing b. This construction corresponds to a two-cutoff 
lattice regulator (TC), where fermion momenta are cut off at a much larger scale 1// than the gauge 
boson momenta 1/6. Doublers are decoupled by introducing a naive Wilson term as in the Roma 
approach p. However, since the chirally breaking effects due to the Wilson term are relevant only at 
scales of the order of the fermion cutoff, gauge boson momenta are not large enough to prove these 
interactions, in contrast with one-cutoff (OC) formulations. As a result, due to the separation of the 
cutoff scales, no hne-tuning is necessary to recover an approximate chiral gauge symmetry, as long as 
the ratio f /b is small enough. In other words, gauge invariance is broken but there is a small parameter, 
f /b, which controls the strength of this breaking. As we will see, there are good arguments to believe 
that a “small” breaking of gauge invariance does not have a dramatic effect at large distances [Q. The 
two-cutoff construction is the hrst proposal that has a control over the strength of this breaking. 

In this paper we want to consider the Wilson-Yukawa picture of the TC formulation of 0 and 
show evidence that in fact a chiral phase exists, in contrast to what was found in OC formulations. In 
order to address these issues numerically, we have considered the simplest relevant model, which is a 
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U{1) compact theory in two dimensions, in the quenched approximation. The model contains one left- 
handed fermion coupled to U{1) gauge helds [|. In order to get rid of doublers, a right-handed fermion 
is introduced, which only couples through the Wilson term. In the global limit {g —>■ 0, i.e. only the 
pure gauge transformations are coupled) the Wilson-Yukawa picture is a scalar-fermion model with 
an exact global symmetry U{1)l x U{1)r. This model has been previously considered in a standard 
OC formulation in [Q, where it was found that the fermion spectrum was either doubled (in the weak 
Yukawa phase) or vector-like (in the strong Yukawa phase). We will show that in the TC formulation 
a truly chiral phase exists. 

The paper is organized as follows. In section 2, we review the relation between gauge-non-invariant 
formulations of chiral gauge theories and models of the Wilson-Yukawa type. In section 3, we present 
the two-cutoff lattice formulation of the two-dimensional model under study and its equivalent Wilson- 
Yukawa picture. In section 4, we discuss the expected phases of this model and why a phase with 
massless chiral fermions should exist. In section 5, we present our numerical results on the fermion 
spectrum, which support the existence of a chiral phase. In section 6 we briefly discuss the effects of 
unquenching and conclude in section 7. 

2 Wilson-Yukawa picture of xGTh 

As explained in the introduction, all proposals to regulate chiral gauge theories which break the gauge 
symmetry are equivalent to a model with an exact gauge invariance and additional charged scalars p[. 
This is simple to see. Let us consider a lattice action that contains gauge-breaking interactions, 

S[U, T] = Sg.ip, T] + SSn.g.p, (1) 

The path integral is 

Z = J VU^ J dTdT » “ (2) 

Since the group is compact we can multiply Z by the volume of the group / VQ, which is an irrelevant 
constant factor: 

Z = j j J dTdT (3) 

Performing a change of variables, from to = Q{xpp)Q^{x + fi) (i.e. the gauge-transformed 
variables under the lattice gauge transformation fl), and from T to = f2(a;)\I'(a;) (only for the 

charged fermions), and using the invariance of the measure and Sg,i. under a gauge transformation]^, 
we get 

Z = 'DUg J dTdT ‘ (4) 

It is easy to see that (^) is gauge-invariant under a new gauge symmetry under which the held 
transforms 

Ug{x) —>• ^{xpp)^"^{x + jl) f2(a;) —f2(a;)<h(x)^ 
p) ^ <|>(a:)T(")(a:) ^ ^(’"^(a:), (5) 

^When the quenched approximation is relaxed more species are needed to cancel gauge anomalies. 

^We are on a lattice, so even if the fermions are chiral the fermion lattice measure is invariant under the unitary 
transformation under Q. 
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where c, n are charged and neutral fermion fields respectively. This is because and are invariant 
under this transformation. This theory is then a gauge theory with extra charged scalars in a non-linear 
realization (hi is unitary). These scalars are the pure gauge transformations, which couple through 
the non-invariant terms in the original action. It is clear that in any successful proposal for regulating 
chiral gauge theories, these scalars must decouple from the light physical spectrum. In the case of a 
spontaneously broken gauge theory these degrees of freedom remain, since they become the longitudinal 
gauge bosons. This is however not our objective here, since we aim at regulating a symmetric chiral 
gauge theory. The interesting case of spontaneous symmetry breaking will be considered elsewhere. 

Foester, Nielsen and Ninomiya (FNN) gave in |P] a good argument why a gauge theory should be 
recovered at long distances in theories like (|^), provided that the strength of the interactions in SSn.g.i. 
is “small”. Their argument is quite simple. Let us suppose that the characteristic coupling of the 
non-invariant terms is arbitrarily small at the cutoff scale. Then in the gauge-invariant picture of the 
theory (^), this implies that the hi fields are very weakly coupled. Consequently, the free energy will be 
maximized when the hi variables are decorrelated at distances of the order of the cutoff (i.e. the lattice 
spacing) or, in other words, these degrees of freedom are effectively very massive. Then it makes sense 
to integrate them out in order to obtain an effective theory at low energies. At distances much larger 
than the cutoff, but still smaller than the correlation length of the gauge-invariant degrees of freedom, 
it can be argued that the effective action should be local (after all it comes from the integration of 
heavy particles) and exactly gauge-invariant (it is clear from eq. (^) that if we perform the integration 
over n what remains is an exactly gauge-invariant theory due to the exact symmetry (|^)). In this 
situation, the only effect of the hi fields is a renormalization of the gauge-invariant coupling in Sg,i,. 

The difficult part in realizing the FNN mechanism is to satisfy the requirement that the interactions 
of the n fields are “small”. The typical situation is that the hi fields are coupled through a higher¬ 
dimensional operator, like a Wilson term, which is suppressed by the fermion cutoff. However, an 
operator of this kind does not satisfy the requirement. This can be seen already at the perturbative 
level. Although the Wilson term is an irrelevant operator, which vanishes naively in the continuum 
limit, it appears in ultraviolet divergent loops, inducing corrections to the scalar action of 0(1) at all 
orders. In contrast, if the boson cutoff (H and U fields) is much smaller than the fermion cutoff, it can 
be shown that the Wilson term is really an irrelevant coupling at long distances, because the gauge 
breaking effects induced by loop corrections are suppressed by powers of the ratio of the boson cutoff 
to the fermion one, to all orders. For details on the power counting arguments specific to the lattice 
TC formulation, we refer the reader to p] and Q]. The FNN conditions are then satisfied by the TC 
formulation if the ratio of these cutoffs is small enough, and we expect that the mechanism is realized 
and an effective gauge invariance is recovered at large distances. 

3 The two-cutoff lattice formulation. 

In 0 and |^, we presented a two-cutoff lattice construction of chiral gauge theories. A different 
proposal in the same direction was given in 0. The advantage of using two-cutoffs as we have 
explained is precisely making the violations of gauge invariance small at large distances. With one 
cutoff a non-perturbative tunning of counterterms would be needed to achieve this ^ . 

The way in which the cutoff separation is implemented is by the use of two lattices. The gauge 
degrees of freedom are the Wilson link variables of an euclidean lattice of spacing 6, Cf,. We will call 
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s the sites of the 6-lattice. The gauge action is the standard one, 

S, = IZEli - + ul.)] 

y s fi<u ^ 

= U^is)U,{s + il)U;:\s + u)U-\s). (6) 

Fermions on the other hand live on the sites of a hner lattice Cf (some integer subdivision of Cb). 
We refer to the /-lattice sites as x. In order to decouple the unavoidable doublers, a Wilson term is 
included in the fermionic action. For each charged chiral fermion, a singlet of the opposite chirality 

is needed. The fermions are coupled to the gauge helds through a standard lattice gauge-fermion 

coupling on the /-lattice. The /-link variables are obtained through a careful interpolation of the real 
dynamical helds, i.e. Uf^{s). As long as the interpolation is smooth, it is clear that the separation of 
scales is achieved in this construction, since the high momentum modes of the gauge helds are cut oh 
at the scale b~^, while the fermions can have momenta of 0{f~^). The lattice action for a charged 
left-handed fermion held is then given by 

+ + - ^{x +il)-ff,[PR + ul{x)PL]'^{x) 

^ o;,// 

+2 d r ^ T(x)\k(a;) — ^ ^['F(x)T(x + fi) +'^{x + /i)T(x)] -|- y 'F(x)T(x). (7) 

X ^ X,fl 

where we have introduced a bare mass term for later use. When we refer to /-lattice quantities we 
use / = 1 units and when we refer to 6-lattice quantities we use 6=1 units for notational simplicity. 
This should create no confusion. 

The TC path integral is then simply 

Z = f VU [ c/TdT (8) 

JCi JCf 

where VU is the standard 6-lattice gauge measure. 

We still have to give a precise dehnition of the u^{x) in terms of the 6-lattice gauge helds h/i(s). For 
a detailed explanation of how to construct such an interpolation for the general case of a non-Abelian 
gauge theory in 4D, the reader is refered to |^. The interpolation is of course gauge covariant under 
gauge transformations on the 6-lattice, i.e. there exists a gauge transformation on the / lattice, uj{x), 
such that 


u^[U^]{x) = u{x) Uf,[U]{x) a;^(x + fi) = u‘^{x). (9) 

In this way, any functional that is gauge-invariant on the /-lattice is automatically gauge-invariant 
under the 6-lattice gauge transformations (notice that the relevant gauge symmetry is that on the 6- 
lattice). Similarly, it is covariant under the remaining 6-lattice symmetries (90° rotations, translations, 
spatial and temporal inversions). This is important to ensure that Lorentz symmetry is recovered in 
the continuum limit 6 —0. Finally, the interpolation has to be as smooth as possible in order to 
achieve the cutoff separation. In , it was shown that a sufficient condition to ensure that the gauge 
symmetry violations induced by the Wilson term are suppressed by //6 is that the interpolated held, 
in the limit / —>^ 0, describes a differentiable continuum gauge held inside each 6-lattice hypercube. 
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and its transverse components are continuous across hypercube boundaries. We refer to this property 
as transverse continuity. 

Now, we want to derive the equivalent Wilson-Yukawa picture. Following the steps of the previous 
section and using the property the model (H) is exactly equivalent to, 

Z = ( vn [ VU ( d-^d^ 

where the a; helds are dehned by (|^) and are coupled uniquely through the Yukawa terms. The oj helds, 
which are determined by the and helds and the interpolation procedure, can be interpreted as 
smooth interpolations of the helds to the /-lattice. In general the transverse continuity property of 
the interpolation also implies that the / —> 0 limit of a; is a continuous held. The high Fourier modes 
q > h~^ of u! are then strongly suppressed. 

Now, in this picture the model has an exact gauge symmetry: 

U^{x) —> ^{x)U^{x)^\x + fi) > f2(a;)<F(a;)'*' 

(j){x)^L{x) (11) 

where 0 is a functional of <F, dehned by the condition, 

u^[U‘^]{x) = 0(x) Uf,[U]{x) (j)\x + ft). (12) 

The new gauge invariance comes at the expense of having unphysical degrees of freedom, uj. Only 
if the u! helds decouple does ([Tl|) imply the true gauge invariance of the original theory (|]). 

3.1 A Two-Dimensional Model 

In order to begin the numerical study of the models (pT]) , we have considered the simplest non-trivial 
model, which is a two-dimensional chiral U{1) theory. The two-cutoh lattice action for the gauge 
non-invariant picture is given by 


S 


^ ^gauge 
s^Cb 


^ ^ ^fermion’! 
x£Cf 


(13) 


with. 


C 


gauge 


c 


fermion 


—^Re Uuis) 
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I E + o:)Pl + (s; + s;)Ph1'|' 


(14) 


where the covariant and normal derivatives are given by Zl+T(x) = M^(a;)\k(a; -f /i) — \k(a;), Zl/\k(a;) = 
4/(a;) — u^fx — {x — p,), d'^ = Zl+|„=i and dp = Dp\u=i. As explained above the link variables 
are interpolations of the real dynamical helds Ufs). The explicit expression is given in the appendix. 
The held Uuis) = Ui{s)U 2 {s + i)17|(s + 2)Ul{s) is the plaquette variable on the 6-lattice. 
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This theory is of course anomalous. In order to get rid of the gauge anomaly more flavours should 
be introduced. However, we will only consider here the quenched approximation and the anomaly is 
not present in this limit. The reason why already the quenched model is interesting is two fold. Firstly, 
in the TC construction, it can be shown that fermion loop corrections, which in general modify the 
scalar interactions, are suppressed by powers of f /b. In this situation, it is expected that unquenching 
is not going to change the phase diagram in any essential way. This is of course not true in the OC 
formulation or in an anomalous theory, where in general fermion loops generate corrections to the 
scalar potential of 0(1). Secondly, Wilson-Yukawa models failed in giving a chiral theory already in 
this approximation, so it is important to check that a different picture emerges with the TC method. 
We will comment on the expected effects of unquenching in an anomaly free model in section 5. 

The Wilson-Yukawa picture is readily obtained from (|iy) and (p!4[) . The action is given by, 

S' = (1 si 

^ j *'^gauge ' / ^ fermion' 

sGCb 


with. 


rwy 

^gauge 


i?e Ui 2 {s) 
9^ 


^fermion 


= 7;j:^'iMD* + D-)Pi, + {a+ + a-)PR\f 

+ + coPl)^ - ^ ( ^u^Pr + ^Pl E d+d;u;^ 

In the case b = f, this model is identical to the one considered in for k = 0. 

As explained before, the model (|l6l) is invariant under a new U{1)l gauge symmetry. 


Ufs) - 

^s) U^{s) <f>\s + fx) 

Ti(x) - 

^ (j){x) Ti(x) 

uj{x) - 

->• Uj{x) (j){xf 


where (j){x) is obtained from ^(s) and from the condition 

ufU‘^]{x) = 4>{x) ufU]{x)(j)\x + fi). 
There is also a global U{1)r symmetry under which 


uj{x) 


(j)R{x) 
(j)R{x) uj{x) 


Finally, a.t y = 0 there is the famous shift symmetry 

'^R —^ d'i? + cr, 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


which ensures that the mass of a fermion with the same quantum numbers as the '^r should vanish 
in the limit y ^ 0. 
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Since we are not interested in a spontaneously broken theory, the a; helds that represent the pure 
gauge degrees of freedom should decouple from the light spectrum. Thus the target theory we would 
like to recover in the continuum limit contains a massless charged left-handed fermion and a free 
right-handed fermion: 

^target + ^{]PPL+ ^Pr)^. ( 21 ) 

As a further simplihcation, we are going to consider the weak coupling limit 0 of this model. 
In this limit the local U{1)l symmetry turns into a global one. Again this has been the approximation 
used in previous studies of standard Wilson-Yukawa models p. In this case, since the Wilson-Yukawa 
model (|T4|) is gauge-invariant, we can simply set = 1 ^ = 1 and what remains is a model of 

fermions coupled to the scalars uj{x). In this case, the uj helds depend only on the fl. The expression is 
given in the appendix, (A.9)~(A.ll). The target theory in the global limit corresponds to a theory with 
free massless fermions. It is easy to trace which of the fermions will couple to the gauge helds, when 
gauge interactions are switched on, by their quantum numbers under the residual global symmetry. 


4 Scenarios for an effective theory 

It is instructive to review what was found in the OC formulation of this model in ref. [Q, in order 
to better understand in what sense the two formulations diher. The model considered in p] is simply 
(H), with an additional kinetic term for the scalar helds: 

2 

^scalar = ^ f2(s + fl) + {s + jl) fl(s) ]. (22) 

^l=l 


The scalar theory (^^, for f = b (OC), is the well-known XY model, which is known to have a 
phase transition at k = Kc ~ 0.56, which separates a vortex phase (VX) for k < Kc with hnite scalar 
correlation length from a spin-wave (SW) phase, where the scalar correlation length is inhnite 


In the quenched approximation, the fermions do not modify the scalar interactions, so this scalar 
dynamics is independent of the scalar-fermion couplings r and y. Since we are interested in a theory 
with no spontaneous symmetry breaking, we should consider the VX phase {k < Kc). Without any 
tuning towards the critical line, the scalars have masses of the order of the cutoff. 

The scalar correlation length decreases as k, is lowered, in agreement with the FNN conjecture, 
since the helds interactions are controlled by k. Alternatively, interpreting helds as the pure 
gauge transformations in a U{1) gauge model, k measures the strength of the breaking of gauge 
invariance at the cutoh scale. According to the FNN argument, at distances much larger than the 
scalar correlation length, the theory should have an ehective gauge symmetry. 

In this model, there are also interactions between the scalars and the fermions through the Wilson 
term in (^). This Wilson term is not small for large r, y and it induces large couplings between the 
fermions and the scalar helds. However, at least naively it seems that the FNN argument can still be 
applied since, in the quenched approximation, the scalar correlation length is controlled only by k and, 
for K -C Kc, it is of the order of the cutoh. At distances that are large compared to this correlation 
length, the physics can be described by an ehective theory where the scalar helds are integrated out. 
The ehective Lagrangian should be local, because it comes from the integration of very heavy particles 
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and it is exactly gauge-invariant, due to the exact symmetry ([T7|). Of course, this is not very useful, 
because it does not tell us what the light fermion spectrum is or if there is any. 

When the Yukawa couplings r and y are small, perturbation theory in these couplings can be 
applied and it is clear that the asymptotic fermion states are the original ones in the Lagrangian. 
Since they have chiral quantum numbers, they should be massless in the FNN regime, already at 
hnite lattice spacing, to ensure gauge invariance. This is precisely what OC studies found in the VX 
(or PM in 4D) weak Yukawa phase, VXW (PMW). In the weak coupling regime, fermion masses are 
proportional to the scalar vev, i.e.. 


m 


{a), 


(23) 


thus vanishing in the symmetric VX phase [k < Kc)- However, as is well-known, the spectrum is 


doubled in this phase, because both doubler and light mode masses behave as (|^) . 

For large Yukawa couplings, the problem is much more difficult. Strong coupling expansions and 
numerical simulations in OC constructions have shown 0, i that in fact there is a free massive 
fermionic spectrum, which gets light when y is taken to zero (the doubler masses however do not go to 
zero in this limit). This is the strong VX (PM) Yukawa phase, VXS (PMS). In this phase, the fermion 
masses are proportional to a globally symmetric condensate of the form. 


m ~ (Re[ f2(s) f2^(s-I-/t)]) 


(24) 


which is non-zero in the VX phase. Although the doublers can easily be decoupled in this phase, the 
light massive fermions were found to transform vectorially under the global group. Obviously they 
have to be different from the ones appearing in the Lagrangian ([T6|), since these have chiral quantum 
numbers. Notice that according to FNN, this is in fact the only possibility: if the light fermions are 
massive at hnite lattice spacing, they must be vector-like in order to satisfy gauge invariance. The 
transition line between the VXS and VXW phases was found to be, for this model at ?/ -|- 2 r ~ 1, 
independently of the individual values of y and r [Q. 

The existence of a strong Yukawa phase implies that, even though the global chiral symmetry is 
not broken (in 2D this has to be the case even in the SW phase according to Mermin-Wagner-Coleman 
theorem), mass terms appear. The solution of this apparent contradiction is that the fundamental 
massive fermions are not the ones in the original Lagrangian (^), but new fermionic composites made 
out of the original ones and the scalar helds, which are expected to form in the strong coupling regime 
It is easy to see that the Dirac fermions. 
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(25) 


transform vectorially under the global group. Notice that for f /b = 1, D = cu. is charged under 
the U{1)l and is charged under the U{1 )r. Since it is conventional to gauge the these 

composites are often called charged and neutral, respectively. It is clear that mass terms for these 
helds do not break the global symmetry. We will call physical the Dirac held made out of the original 
helds in (^, i.e. = 4/^-1- since it is the one we expect to hnd in the light spectrum according 
to (0). 

When the couplings y and r are large (0(1)), the fermion-scalar couplings in (pffil) are strong enough 
to favour the formation of the mirror composites, uj and The numerical simulations of |^, ^ 
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showed that the strong phase could be reached by making r large while keeping y arbitrarily small, 
as long as the condition y + 2 r > 0(1) was satisfied. In fact, the numerical results for the neutral 
propagator (d'O)) were nicely reproduced by the first-order corrections in the hopping parameter 
expansion 0], which does not distinguish between the regions r ^ y = 0(1) and y <^r = 0(1). 

To first order in this expansion, the neutral field is a free Wilson fermion. The Lagrangian for the 
neutral field is given by, 

r - 2 

+ y f ^ (26) 

2 ^=1 

The zeroth-order propagator is given by. 




-i Ea. 7m sin q^,PR + M{q) 

M{qy 


(27) 


with M{q) = 1 — a r cos q^ and a = l/{y + 2r). At higher orders, the propagator can be writen in 
general as 


S^^\q) = 


i sin q^{f^^\q)PR + f't\q)PL) + M^'^){q) + M^^\q)-i^ 


p(0 


O)/ 


(28) 


At first order in one gets |l[] 

fniq) = 1 fL\q) = 

M^^\q) = M{q) M^^\q) = 0 

where 2^ is the well-known condensate 


(29) 


2^ = (Re([a(n)[a^(n H-/i)))n- 


(30) 


In the strong Yukawa phase, plays the same role as the vev in the broken phase. However, it is 
invariant under the global chiral symmetry, as opposed to the normal vev. This is of course the reason 
why it induces fermion masses in the symmetric phase. It is easy to see that the propagator defined 
by (1^) and (^9]) is equivalent to a free Wilson propagator with. 


y 

Veff = '^eff = -1 

/2 


ZlR — 2 : Zri — 1, Zll — Zrr — z ^. 

The Lagrangian for the charged fermion is given by, 

r - 2 

2 ^=1 


(31) 


(32) 
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In this case, the numerical results were not well reproduced by the hopping expansion. On the contrary, 
there was strong evidence that this fermion is not a free asymptotic state but most probably a 
two-particle state, 

In summary, OC studies found that both the strong and weak Yukawa phases failed in giving a 
true chiral spectrum. As we have seen, these results are in nice agreement with the FNN mechanism. 
The scalar helds decouple at large distances and an effective gauge-invariant theory remains. However 
the gauge symmetry can be realized either through a massive vector-like fermion spectrum or through 
a doubled massless one. FNN does not exclude the possibility that composite fermions exist at large 
distances nor that there is doubling. We will now argue that in the two-cutoff formulation (pTl), there 
is a richer structure in the strong Yukawa phase, depending on the individual values of y and r. 

Clearly the only difference between OC and TC is in the scalar dynamics. In the TC case we have 
set K = 0 for simplicity, so we are deep in the VX phase. In OC studies the fermion spectrum has 
not been analysed in detail in this extreme region of phase space. Probably one of the reasons is that 
it is much more demanding numerically. The other reason is that in the unquenched theory fermion 
loops introduce large (typically 0(1)) corrections to k, so that in the full theory it makes no sense to 
consider the case k = 0 unless some hne-tuning is done. In the TC case, as we will see, fermion loops 
introduce only a small correction to k, provided there are no gauge anomalies. 

The main difference between the two formulations is the fact that, due to the interpolation, the 
high Fourier modes of the scalar held are suppressed. Let us consider a region of phase space with y 
arbitrarily small and r = 0(1). We use / = 1 units unless stated otherwise. According to the OC 
studies this would correspond to the strong phase, where the spectrum is vector-like. However, for these 
couplings, the light lattice mode (i.e. p —0) is coupled to the scalar held only through a perturbatively 
small y Yukawa coupling, and a higher-dimensional operator. At small lattice momentum, the ehect 
of this operator is small. Notice that this is true only for TC. In the OC case, the Wilson term also 
induces a large 0 {r) coupling between the light mode and the doublers, so it is not possible to treat it in 
perturbation theory. In the TC case, this coupling is small, because it involves high-momentum modes 
of the scalar held (~ /~^), which are suppressed. On the other hand, for the doubler poles, p —> tt, 
the Wilson term is ehectively a Yukawa coupling of strength 0(r), in both OC and TC formulations. 
In this situation, we expect that at small enough y, the light mode enters into a weak phase, because 
the strength of the Yukawa couplings is small at p —> 0, and there is no reason to believe that the 
formation of the mirror state is dynamically favoured. On the other hand, the doubler modes could 
remain in the strong phase, since for them the Yukawa couplings are of 0(1). This of course would be 
a most desirable scenario, since the doublers would be decoupled from the light spectrum, while the 
light modes would remain massless without any extra tuning, in the same way as they are in the weak 
Yukawa phase. In other words, the light fermion spectrum would be chiral and not doubled. For this 
reason we will refer to this phase as the chiral phase. Clearly, the fact that the scalar held momenta 
cannot get as large as the fermion cutoff (the doubler momentum) is essential in this argument. 

If a new phase of the kind discussed above exists at small p, it cannot be described by the low¬ 
est orders of the hopping parameter expansion and it is reasonable to expect that this expansion 
should signal its failure to converge at small y. The hrst-order hopping result (^) corresponds to 
the resummation of a set of one-loop tadpole diagrams, which is a geometric series in the quantity 
(az)"^ S/j(sinq'^)^/M(g)^. If the corrections that are neglected at the next order are of the same order 
of magnitude as the second term in this series, the hopping expansion is expected to be a good approx¬ 
imation only if (az)^ X]^i(sin g^)^/M(g)^ -C 1, which is a momentum-dependent condition. At small 
momentum, M(g) —1 — 2ar ~ y/2r and the convergence condition would require zj{yL) -C 1, where 
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L is the lattice size. For small enough y the condition is clearly not satished, since is independent 
of y and, as we will see, in the TC case, it is ~ 1. 

In the region of small y and at small momenta, we expect to hnd a chiral phase, in which the 
light fundamental helds are those in the original Lagrangian (pIBf) that we have called physical. This 


phase is clearly more difficult to analyse than both the strong and weak phases, since doublers and 
light modes behave in a very different way, even though they are excited by the same held. At large 
momenta, the hopping expansion is expected to be a good approximation, while at small momenta 
weak perturbation theory should be a good approximation. The exact result for the propagator of the 
physical helds in momentum space is then expected to be an interpolation between the hopping result 
at large momenta and the weak expansion result at small momenta. 

At small momentum and to zeroth order in y and r: 


(SW(«))° = 


>11 r^sinq^) 


-1 


Q 


0 . 


(33) 


Higher-order corrections can induce a renormalization of the helds so in general we expect that the 
propagators will be well reproduced by 


S''’>(9)ih = = 


* E/9 7m sin Q/q ‘ 






i E/q 7m sin q^q 


% 


0 . 


(34) 


The LL, RR components of are zero, because there is no chiral symmetry breaking. However 
all the remaining components of the neutral and charged propagators, i.e. and 

are non-zero. In weak perturbation theory, they correspond to multiparticle states made of various 
combinations of the massless helds, 4/x^ and the scalar uj. For instance, at large times, we expect 

(n) fc) 

that the propagation of and 4/^ is dominated by two-particle states, i.e. 

SLRix, y) ~ y) Guq{x, y), y) ~ y) y), (35) 


where is the scalar propagator. Since these propagators involve scalars, which have masses of the 
order of the boson cutoh, they should decouple from the light spectrum. 

At large momentum, a better approximation is provided by the strong coupling expansion. In 
this case the physical propagators S'jfx and combine with the remaining neutral and charged 
components to form massive Dirac helds. Thus we do not expect light doubler poles. The physical 
propagators at large momentum should behave as 


RL,LR 


Zr,l 


-i Em 7m sin + Mr^l 
EM(singM)^ + ^l,L 


% 


TT. 


(36) 


To summarize, the expected phase diagram in the symmetric phase, in the TC formulation, has 
three regions (not necessarily separated by phase boundaries). At small ?/, r -C 1, there is a weak 
phase where all fermions are massless, doublers and light lattice modes. The ehective theory in this 
regime is diherent from the target theory (^) in that there is a doubling: for each mode with a Pl 
coupling to the gauge held, there is one with a Pr coupling. At large y,r = 0(1), there is a strong 
phase where all fermions are massive, although chiral symmetry is not broken. The real asymptotic 
fermion state is a singlet fermion (the charged one is possibly just a two-particle state, as favoured 
by the OC studies [^), which only couples to the gauge helds through higher-dimensional operators. 
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This is also far from (pT]). Even if the charged fermion existed as a free fermion in this phase, it would 
couple vectorially to the gauge fields, so again the continuum effective theory would be different from 
the target theory. Fortunately, there is a third region at y -C r = 0(1), which was not present in OC 
models, where the doublers are decoupled from the light mode, but the latter is only weakly coupled to 
the scalar field and remains in the weak phase. This chiral phase is expected to give rise to the target 
theory defined by (^). Our main purpose in the following is to show clear evidence of the existence 
of this phase in the TC formulation. 


5 Numerical results 

We have computed the fermion propagator for the different fields and by inverting the 

appropriate fermion matrices on configurations of the uj fields that were obtained by interpolating a 
set of completely random f2(s) configurations, in the way described in the appendix. This corresponds 
of course to the case k = 0. For the fermion fields, we have used periodic boundary conditions in the 
spatial direction and antiperiodic in the temporal direction. This is necessary to study a regime in 
which we expect to find massless fermions. On the other hand, the scalar fields O are periodic in both 
directions. 

With two lattices, there is some freedom in the choice of the ratios b/f and L/h and it is not 
obvious how the continuum limit should be defined in this case. Clearly, the two cutoffs must be sent 
to inhnity, i.e. f,b^0. When gauge interactions are switched on, we expect to have a mass gap (.^), 
so that the continuum limit is defined by. 


b 

1 


0 


(37) 


This however might require a simultaneous tunning of //6 —0 to ensure that we do not exit the FNN 
regime. In |^, it was argued that a conservative way of taking the continuum limit without exiting 
the FNN regime (i.e. ensuring that fermion loops contributions to the scalar dynamics are suppressed) 
would be 


L k 

b b 


0 , 



(38) 


however, it is not possible from theoretical considerations to determine quantitatively how these two 
limits should be correlated. To understand this issue a full simulation with gauge interactions would 
be required since, in the simplified model we are considering here, there is no physical mass gap. 

We have used units of / = 1 throughout the paper, unless stated otherwise, for notational conve¬ 
nience. However, it is important to remember that the boson cutoff is 1/6 and thus in /-lattice units, 
masses of 0(//6) are of the order of the gauge boson and scalar cutoffs. Finite-volume effects can be 
studied by keeping the couplings fixed in 6 units, while increasing the ratio L/b. This might require a 
simultaneous change in the ratio 6// to ensure that we remain in the FNN regime. We will refer to 
the different lattice sizes with the notation {L/b)b/f. For convenience we have used square lattices, so 
that L/f is the same in both space and time directions. 
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Figure 1: Scalar Propagator in 
ansatz (|n) at small momenta, 
polynomial in 


momentum space in a Ss lattice. The dash line is the ht to the free 
The solid line corresponds to a £t to the inverse of a fourth order 


5.1 Scalar Dynamics 

The main difference with OC simulations is clearly in the scalar dynamics. The scalar propagator 


G.{p) 





(39) 


is quite different from the one corresponding to a free massive boson held. In Fig. 1 we show the scalar 
propagator in momentum space for an 83 lattice. The two sets correspond to Guj{qx,0) and Gi^{0,qy) 
respectively. They are compatible within statistical errors, as expected from the lattice 90° rotational 
symmetry. At small momentum the propagator can be htted to the free ansatz (dashed line in Fig. 
1 ): 


GM) 


q"^ + ml' 


(40) 


with = 2 (1 — cosg). In this case, the htted scalar mass is roughly = 1.7 in b units. This is 
in agreement with the naive expectation that this mass is related to the boson cutoh 1/6. In fact, 
if the interpolation was strictly local (i.e. the interpolated held in each b plaquette depends only on 
the surrounding links), the correlation length of the ut helds at distances larger than 2 b would vanish. 
In this case, it would not be appropiate to talk about a scalar mass. However, due to the fact that 
the group is compact, a smooth interpolation cannot be strictly local 0. Empirically we hnd that 
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Figure 2: Scalar propagator in the large-momentum region for two values of//6=l/8 (triangles), 1/4 
(circles). The curves are fits to the inverse of a fourth-order polynomial in in the whole momentum 
range. 


the non-locality of the interpolation is nevertheless consistent with an exponential decay of the scalar 
propagator at large distances (> 2 b), with a mass of the order of the boson cutoff. 

At large momentum, however, the a; propagator differs considerably from the free ansatz. Since 
the interpolation is smooth, the momentum modes higher than g ~ //& (boson cutoff) are expected 
to be suppressed for small f /b. Empirically, the inverse propagator at large momenta can be htted 
for these lattice sizes to a polynomial in g^ of fourth order (solid line in Fig. 1). In Fig. H, we show 
the large-momentum behaviour of the scalar propagator for two different values of the two-cutoff ratio 
f/b. Clearly, the high-momentum suppression increases as this ratio is decreased, as expected. 

Concerning the condensates, there is also some difference with the OC case. The chirally breaking 
condensate, 

^ = (41) 


is expected to be zero, since we are in the symmetric phase. Numerically, it is zero within statistical 
errors. 

The condensate relevant for the strong phase. 




2 


'2L2 


y/ Re[ uj{x)u\x -|- /i) ]), 

X,fl 


(42) 


is plotted in Fig. 3 for two lattice sizes of L/b = 8, 4, as a function of the ratio f/b. It is easy to show 
from the explicit expression of the interpolated uj helds that. 


lim z‘^ = 1 (43) 

f/b^O 
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Figure 3: condensate for two lattice sizes L/h = 8,4 as a function of the ratio f /h. The lines are 

quadratic hts. 


//& 

L/h 

{h units) 

1/8 

6 

1.70(6) 


8 

1.72(6) 


10 

1.75(15) 

1/4 

8 

1.69(6) 


10 

1.54(10) 


Table 1: Scalar mass in units of b for various lattice sizes from momentum hts (PUI). 


for hxed L/h. This is in contrast to what happens in the OC formulation, where = 0 for k = 0. This 
important difference between the OC and TC cases is not strange, since gets important contributions 
from short distances (as opposed to n, which only depends on the long-distance scalar dynamics) and 
the short-distance scalar dynamics is drastically modihed by the interpolation. 

One important question is what the volume dependence of these results is. Since the O helds are 
random at the scale b we do not expect an important dependence on L/bf^ We have computed the 
scalar propagators for various lattice sizes, and we hnd that the scalar “mass” (as obtained from the 
momentum £t (|4DD) is proportional to the ratio f/b, as expected. This means that the scalar mass 
in b units does not depend on f/b. On the other hand, there is no clear dependence on L/b. This is 
shown in Table 1. In all the cases the scalar masses are of 0(1) in 6-lattice units. 

The volume dependence of the condensate z^ is shown in Fig 3. For hxed f /h, there is some 

^Unfortunately this cannot be made into a rigorous statement because of the winding fields N{s), which can grow 
with L/h. See the appendix. 
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dependence on the volume^. In order to keep constant, the inhnite volume limit must be taken in 
a correlated way with the ratio j jh. A sensible way of taking the inhnite volume limit might be by 
following lines of constant z^ in the {f/b,b/L) plane. From Fig. 3 we see that, according to this rule, 
f /b does not need to decay as fast as b/L. 

5.2 Fermion Spectrum 

The fermion propagators have been obtained with the conjugate-gradient method for several lattice 
sizes and y couplings. Since we are only interested in the regime where doublers are decoupled, we 
have kept the Wilson coupling hxed at r = 1. According to the OC phase boundary y + 2r = 0(1), we 
are then in the strong phase. As we will see, the strong phase in the TC has more structure, depending 
on the value of y. 

The various propagators are given by 

St,{t) = {Y. >n. (44) 

xi,yi 

where t = 1,..., L, the index a refers to the neutral (n), charged (c) or physical (p) fermion for pi = 0, 
and the corresponding spatial doublers (nd), (cd) and (pd) for pi = tt. i, j = R, L for the different 
chiralities. For every inversion the time slice at the origin X 2 is chosen randomly. The number of 
sampled scalar conhgurations is typically of 0(2 — 5 10^). 

As in the OC construction, invariance of the action (|26|) under PT, i.e. a spatial and temporal 
inversion^, implies 


^“(f) = -75 5“(T-t) 75, (45) 

which means that there is antisymmetry under t ^ T — t for the chirality-breaking components, RR 
and LL, and symmetry for the chirality-preserving ones, RL and LR. Notice that there is no symmetry 
under T or P, since the action explicitly breaks parity. 

In order to check that these symmetries are well satished we have monitored the behaviour t <T/2 
and t > T/2 of the propagators. The expected PT symmetry is well satished within statistical errors. 
If the low lying spectrum consists of free massive fermions, P and T symmetries are also expected to 
be satished ehectively, while if the spectrum is chiral P and T are broken explicitely. 

5.2.1 Large-?/ Regime 

At large y S> //5, we hnd that all the fermions are massive, as was also found in OC studies. Of 
course this phase has no physical interest since the fermionic states are more massive, m ~ ?/, than 
the scalars and this implies that there is no FNN regime to look at. All the particles, scalars and 
fermions, have masses of the order of the cutoh. In any case, it is interesting to look at this regime 
in order to compare with OC studies and to conhrm that also in the two-cutoh formulation, for large 
Yukawa couplings, fermions are massive even if chiral symmetry is not broken. Figure ^ shows the 
inverse neutral propagators in momentum space for ?/ = 1.5 in an Sg lattice. The solid lines correspond 

^This effect is also due to the winding fields N{s). 

^Notice that T is not equivalent to time reversal in the continuum. 
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Figure 4: ^ at y = 1.5 and a lattice size of 83 . The solid lines correspond to the hts to the 

free Wilson fermion formulae of eq. (4B). The parameters of the hts are given in Table 2. 


to the hts to the free Wilson formula (for zero spatial momentum), i.e. 


SiR,RL(«) 
d”L’.HR(9)'‘ = 


sing 

(sing)2 + + - cosg))^ 

Veff+ ref f {I-cos q) 

(sing )2 + + reff {I - cosg))2 


(46) 


with Zll,rr = All propagators can be htted with a very good to the free Wilson 

propagator. The hopping expansion is a very good approximation at large momentum (g > vr/fo) as 
can be seen in Table 2. We present the results for ge// and rg// of the hts to (]4^), excluding the points 
near zero momentum, together with the result of the hopping expansion to hrst order (|^) . If the 
points near zero momentum are included in the ht the increases considerably, indicating that the 
free propagator (^) is not a good approximation in the low-momentum region. 

This ehect is more clearly seen in the behaviour of the propagators at large times. In Fig. |^, we 
present the LR component of the neutral propagator for two large values of y together with the result 
from the momentum hts. Clearly, there is a systematic diherence at large times, which only decreases 
slowly as y increases. The time propagators can be very well approximated by a two-exponential ht, 
as shown in Fig. ^ but we do not have a deep understanding of the origin of the lighter states that 
seem to propagate at large times in this channel. What is clear from our data, however, is that at 
small times, typically t < b — 2b, the propagators are very well reproduced by the hopping result and 
are thus dominated by the propagation of a free massive fermion. In the future we will explore the 
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y 

Prop 


RR 

RL 

LR 

LL 

Hopp. 

1. 


z 

Veff 

reff 

1.084(5) 

1.090(6) 

1.081(6) 

1.006(23) 

1.081(18) 

1.090(11) 

1.175(19) 

1.087(11) 

1.083(8) 

1.086(5) 

1.093(6) 

1.082(5) 

1.085 

1.085 


z 

Veff 

^eff 

1.084(5) 

1.085(5) 

1.084(5) 

0.992(26) 

1.064(14) 

1.086(14) 

1.176(10) 

1.085(5) 

1.084(5) 

1.085(5) 

1.086(5) 

1.085(5) 

1.085 

1.085 

1.2 


z 

Veff 

reff 

1.082(5) 

1.304(5) 

1.080(5) 

0.997(23) 

1.289(22) 

1.085(14) 

1.172(11) 

1.301(6) 

1.082(5) 

1.085(5) 

1.307(6) 

1.082(5) 

1.302 

1.085 

l^(nd) 

z 

Veff 

reff 

1.084(5) 

1.301(6) 

1.084(5) 

0.969(32) 

1.266(20) 

1.072(20) 

1.175(10) 

1.301(6) 

1.084(5) 

1.085(5) 

1.302(6) 

1.085(5) 

1.302 

1.085 

1.5 


z 

Veff 

reff 

H 

1.009(15) 

1.624(13) 

1.093(7) 

1.175(8) 

1.628(5) 

1.083(4) 

1.085(3) 

1.631(5) 

1.083(3) 

1.627 

1.085 


z 

Veff 

^eff 

1.085(3) 

1.628(5) 

1.085(3) 

0.991(25) 

1.608(19) 

1.083(13) 

1.178(13) 

1.629(10) 

1.085(6) 

1.085(3) 

1.628(5) 

1.085(3) 

1.627 

1.085 


Table 2: Fits of the neutral propagators (light and doubler) in momentum space to the free Wilson 
formula for large y in an Sg lattice {z~^ = 1.085). 
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Figure 5: propagators in time for y = 1, 1.5 and a Sg lattice. The lines following the data points 

are the result of a two exponential £t. The dotted and solid lines, which deviate at large times, are 
the results of the momentum hts of Table 2. 

possibility that the deviations from the hopping result could be explained by the presence of fermion- 
scalar interactions. As explained in ||^ an effective Lagrangian for the neutral held compatible with 
the global symmetry is 

Ceff ^ + mn)¥^^ + A - ¥d,uj). (47) 

OC studies did not find evidence for such interactions. In fact, when the fermions are much lighter 
than the scalar helds, this is to be expected from FNN arguments, since the scalars should decouple 
from the light spectrum, but when the fermions are more massive than the scalars as in our case, there 
is no reason to believe that these interactions are absent. 

The results for the neutral doubler are also shown in Table 2. As expected, there is good agreement 
with the hopping result in momentum space. The large-time behaviour of these propagators, however, 
also shows an anomalous behaviour: the effective mass at large times clearly decreases to a smaller 
value, see Fig. ^ In the case of the doublers, a scalar-fermion interaction of the form (|47|) means that 
the doubler is not stable. It can decay into a light mode and two scalars, since it is heavier. This 
coupling is quite weak because in such a decay the scalars carry large momenta and, as we have seen, 
the large Fourier modes of the scalar helds are strongly suppressed. Empirically, the overlap with the 
lighter states is very small: Z ~ 10“®. 

To summarize, we can conclude from our data that the hrst order hopping expansion provides 
a good approximation for the neutral propagators at large momenta and small times, in the large-?/ 
region. At large times, the ehective mass decreases and it is not clear whether there is an asymptotic 
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Figure 6: propagators in time for y = 1, 1.5 in a Sg lattice. The dotted and dashed lines, which 

deviate at large times, are the results of the momentum hts of Table 2. 

free fermion in this channel. Notice that the FNN argument is not at work here, because the fermion 
mass is much larger than the scalar correlation length. In any case, the fact that the propagators 
are well reproduced by the hopping expansion at small times {t < b — 2b) implies that, by choosing 
f/b small enough, the correlations can get negligibly small at t ~ fe and are not expected to affect 
long-distance observables, if a light sector was present in the theory. 

The charged propagators, both doubler and non-doublers, also behave as massive Dirac fermions. 
However, in this case the hopping approximation is not a good approximation even at small times. 
This is similar to what was found in OC studies. In Fig. we show the RL components of the charged 
propagator, together with a two-exponential £t. In Table 3, we compare the effective masses obtained 
from a single-exponential £t in the large-time region of the neutral and charged LR propagators. The 
difference between them is of the order of the scalar mass from Table 1. This is consistent with the 


y 

(/-units) 

(/-units) 

1. 

0.538(4) 

0.745(27) 

1.2 

0.653(7) 

0.868(43) 

1.5 

0.816(3) 

0.982(21) 


Table 3: Effective masses from a single-exponential £t at large times of the neutral and charged LR 
propagators for large ?/, for a lattice size of Sg. 
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Figure 7: propagators in time for y = 1, 1.5. The lines following the data points are the result of 

the two-exponential fits. 

charged propagator being a two-particle state, as favoured by OC studies. 

5.2.2 Small-?/ Regime 

As we decrease ?/ < 1, the different chiral components of both the light neutral and charged propagators 
start to differ, signalling a violation of parity. Two components start getting lighter than the others: 
^RL = ^RL (i.e. the expected physical helds). From our data we cannot establish 

a transition line in which parity violation shows up. It seems that the behaviour of the propagators 
changes in a smooth way, the violations of parity being present also at large y and getting larger 
smoothly as y is decreased. 

There is, however, a critical value of yc below which the two lighter components and get 
actually massless at finite lattice spacing. This is the onset of the announced chiral phase. It is not 
clear what determines the value of yc- There is one obvious candidate, which is the boson cutoff, i.e. 
f /b in /-lattice units. Although our numerical results are not dehnite on this point, they seem to 
indicate that yc is related to this scale. The physical picture behind this expectation is nothing but 
the FNN conjecture that at distances larger than h/f in /-units, the scalars should decouple. The 
correlation length of the light fermions in the strong phase is ~ ?/, so when y < f /b the scalars decouple 
from the light fermions and a chiral phase appears. 

As explained in the previous section, in order to study this phase, it is better to go to momentum 
space, where the limits g —> 0, tt of the propagator are known. In Fig. we present the RL and LR 
components of the neutral and charged inverse propagators for y = 0.05. Clearly the and 
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Figure 8: Inverse propagators (5'^”^’^"''^^) ^ and (5'^'^^’^'^'^)) ^ in momentum space in an Sg lattice for 
y = 0.05. 


Slr = propagators have poles only at g —> 0, as expected from massless undoubled fermions. 
On the other hand, the components and show no poles, in agreement with them being two 
particle states composed of a massive scalar particle and a massless fermion 

In Fig. 1^, we present a zoom of the massless components in the small momentum region for several 
values of y. The continuum curves are hts of the form (the spatial momentum being zero). 


q{n,c) 

‘^RL,LR 


(?) 


Zr,l ^ sin(g) 

sin(g) '^’^sin(g)2 + 


(48) 


The second term corresponds to higher states, which should decouple in the continuum limit. Their 
masses are expected to be of the order of the boson cutoff or larger. It is clear from the hgures that 
the contamination of the charged chiral mode is larger than that of the neutral one. We show the 
results of these hts for several lattice sizes in Table 4. The volume effects as L/b changes, for hxed //6, 
are quite small and it seems clear that the contaminating mass m' does not decrease as the volume is 
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Figure 9: ^ and (S^lr) ^ 


at small momentum for several values of y in the chiral phase. 


Size 

Zr 

Zl 

m'j^ (&-units) 

^'l 

Zr. 

66 

0.994(1) 

0.543(12) 

1.16(17) 

0.384(25) 

68 

0.994(1) 

0.668(17) 

1.27(34) 

0.257(23) 

84 

0.994(3) 

0.168(33) 

0.95(22) 

2.245(54) 

88 

0.996(1) 

0.494(21) 

1.08(13) 

0.603(84) 

108 

0.995(1) 

0.485(40) 

1.22(44) 

0.52(12) 


Table 4: Z'j^/Z^ and for various lattice sizes, y is kept hxed to 0.2 in units of h. 
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Size 

{b units) 

{b units) 

66 

2.26(11) 

2.08(78) 

68 

2.19(14) 

2.03(11) 

§8 

2.46(13) 

2.32(10) 

lOs 

2.59(51) 

2.03(19) 


Table 5: and from a 

kept fixed to 0.2 in b nnits. 


fit to for varions lattice sizes, in nnits of the boson cntoff b y is 


increased. On the other hand, for hxed L/6, the contamination increases considerably for the charged 
propagator, when the cntoff ratio is increased. 

In Fig. 0 we zoom the low-momentnm behavionr of and It is important to make sure 
that these propagators do not develop poles as y gets smaller or as the continuum limit is approached, 
rendering the theory vector-like. It is not clear how to £t these propagators at zero momentum, since 
they are expected to be multiparticle states. We have used the simplest ansatz. 


n{n,c) 

^LR,RL 


(?) 


sin(g) 


‘■* sm ((/)2 + 


(49) 


which empirically gives a good £t near g = 0. In Table 5, we present the masses obtained from these 
hts for various lattice sizes. It is clear that these masses are of the order of the boson cutoff, as 
expected from (P^D. Also, there is no sign of the masses getting smaller as the volume increases or as 
y is decreased. 

For the neutral doublers, we hnd that the large-momentum (g > vr/fe) and small-time {t < b) 
behaviour corresponds to that of a free massive fermion, which supports the picture that, due to the 
strong coupling between the doubler modes and the scalars, the neutral composite fermion (^5[) is 
formed. In Table 6, we present the hts of the neutral doubler inverse propagators to the free Wilson 
formulae for large momentum (g > n/b). The effective r is in perfect agreement with the hopping 
result. The effective y on the other hand differs (especially for the RL component). This could be due 
to a low-momentum bias, since the hopping result is not a good approximation for momenta of the 
order of the boson cutoff. 

At large times or small momenta, there is evidence of contamination from lighter states, which are 
probably multiparticle states. This is not so surprising since the neutral composite is not expected to 
be stable in this phase. Possible effective interactions, compatible with the symmetries that connect 
the light and doubler sectors are given by 


Qn { ¥P'>uj^Pii -1- dliuj}. (50) 

In the presence of these interactions, the neutral doubler channel at large times will be dominated 

by two-particles states. The decay of the neutral doublers into lighter particles however involves high 
momentum modes of the scalar held, which are suppressed for TC, so we expect that they are relatively 
long-lived and thus dominate the propagators at small times. Notice that in the OC case, this does 
not have to be the case. 


24 








o 


% 



qt 


Figure 10: ^ and ^ at small momentum for several values of ?/in the chiral phase. The lines 

are fits to (H), which give = 2.46(13), 2.46(35), 2.10(55) and = 2.32(10), 2.19(21), 2.07(29)., 
in increasing order of y. 


y 

Prop 


RR 

RL 

LR 

LL 

Hopp. 

0.025 

g{nd) 

z 

Veff 

^eff 

1.083(6) 

0.033(5) 

1.082(5) 

0.992(16) 

0.004(4) 

1.090(8) 

1.173(8) 

0.030(4) 

1.082(4) 

1.087(3) 

0.033(1) 

1.086(3) 

0.027 

1.085 

0.05 


z 

Veff 

^eff 

1.084(5) 

0.061(3) 

1.082(4) 

0.993(16) 

0.024(4) 

1.090(8) 

1.173(8) 

0.057(1) 

1.082(4) 

1.087(2) 

0.059(5) 

1.086(3) 

0.054 

1.085 

0.08 


z 

Veff 

^eff 

1.082(4) 

0.092(1) 

1.081(4) 

0.997(17) 

0.056(5) 

1.092(8) 

1.172(8) 

0.089(1) 

1.082(4) 

1.087(4) 

0.092(1) 

1.086(4) 

0.087 

1.085 


Table 6: Fits of the neutral doubler propagators in momentum space (g > 7r/6) to the free Wilson 
formula for small y in an 83 lattice. 
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Figure 11: Effective mass of the neutral (up) and charged (down) doublers for several lattice sizes and 
y = 0.025. 


In Fig. 1^ we present the effective masses in the RL and LR neutral and charged doubler channels 
for t < b — 2b, for three values oi L/b and hxed f /b. The neutral effective mass has a nice plateau 
at small times. This is consistent with the assumption that this channel is dominated by the neutral 
composite of (|^) . On the other hand, the charged doubler does not show a plateau. The effective 
mass is of the order of the fermion cutoff at small times, but it is nowhere constant. This indicates that 
the charged doubler composite is not formed and multiparticle states are dominating this channel. 

Although our statistics is not good enough to attempt a quantitative description of the doubler 
propagators at large times, the fact that, at small times t < b — 2b, the propagators decay as particles 
with masses of the order of the fermion cutoff implies that, after a few /-time slices, the correlation in 
the doubler channels is too small to influence the long-distance physics. In the simplihed model we are 
considering here there is no more long-distance physics than the free propagation of massless fermions. 
In Figs. 0 and we compare the decay in time of the various propagators, neutral and charged 
respectively, for several lattice sizes. The massless components have constant correlation as expected. 
All the remaining propagators are not simple at large times, but they are clearly split from the light 
sector. For f /b = 1/8, they are suppressed by at least 2-3 orders of magnitude. We notice that the 
effect of increasing f/b is quite dramatic for the charged propagator (see last hgure in Fig. |^. In 
general the splitting between the doubler and non-doubler sectors is larger for the neutral channel. 

To summarize, we have found that at small enough y there is a chiral phase, in which the light 
spectrum consists of two massless fermions with chiral quantum numbers. Their propagators in mo- 
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Figure 12: Number of conjugate-gradient iterations in the inversion of the neutral propagator as a 
function of y in an Sg lattice. 

mentum space show only one pole at g —0, as expected from undoubled massless particles. The 
neutral doubler channel seems to be dominated by a massive fermion at small times, while at large 
times it is probably dominated by multiparticle states. The charged doubler channel seems to be 
everywhere dominated by multiparticle states. These states contain scalars, whose correlation length 
is of the order of the boson cutoff 6“^; thus they are expected to decouple in the continuum limit. 
In order to show that this is indeed the case, it will be necessary to consider the case where gauge 
interactions are switched on, since only in this case is there a physical mass gap that can be used to 
define the continuum limit of the lattice model. We have not found important volume effects as L/b 
is increased, although there is a clear dependence on f/b: the splitting between the light sector and 
the heavy one becomes smaller as f /b is increased. This is a clear indication of the importance of the 
cutoff separation. 

We have not analysed in detail the intermediate region of couplings yc < y < 0(1), and in particular, 
we have not addressed the important issue of the nature of the transition at yc- More statistics will 
be needed for this. However, we have monitored one observable that has proved useful in the location 
of crossover or transition lines in previous studies of Wilson-Yukawa models |]I|. This is simply the 
number of conjugate-gradient iterations (NCG) needed in the propagator inversion for a hxed precision. 
The naive expectation is that the required number of iterations should increase as y decreases, because 
the lowest eigenvalue of the fermion matrix decreases with y. This changes however when we enter the 
chiral phase, since the light mode is massless and the lowest eigenvalue is fixed by the IR cutoff, which 
is the lowest lattice momentum (remember we use AP boundary conditions in the time direction). 
Thus the number of iterations should not keep increasing for y < yc- In Fig. 13. we show the NCG 
for different values of y. There seems to be a maximum very close to the value where we believe yc 
should be, according to the behaviour of the propagators. The position of the maximum is in fact near 
y = 0.15 in / units, which is not far from f/b = 1/8, as naively expected. 
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Figure 13: S^l (full circles), (triangles), (circles), (squares) for several lattice sizes: Gg, 
Sg, lOg and 84, for y = 0.025. 
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Figure 14: (full circles), (triangles), (circles), (squares) for several lattice sizes: Gg, 

Sg, lOg and 84, for y = 0.025. 
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6 Effects of unquenching 

It is clear that unquenching in the OC case can have dramatic effects. By naive power counting, 
fermion loops can induce 0(1) corrections to the potential of the lo fields. In order to be in the FNN 
regime, tunings of k and probably other relevant operators of dimension 4 that are generated at higher 
orders might be needed. 

As explained in |^, in the TC case the situation is very different. Due to the suppression of the 
high-momentum modes of the scalar helds, it can be shown that the corrections induced by fermion 
loops on the scalar potential are down by powers of the two-cutoff ratio|^, 0(//6)^. For this reason it 
is expected that unquenching is not going to change the phase diagram of the fermion-scalar model in 
any important way. A way to quantify this is by looking at how the phase of the fermion determinant 
for the action (pA]) changes under a gauge transformation. The magnitude of this change is a measure 
of the strength of the couplings of the scalar helds a;, induced by the fermion loops. At this stage, 
however, we have to enlarge the fermion content of the model we are considering in order to cancel 
gauge anomalies. Otherwise the changes in the phase are of 0(1) and the two-cutoff method is of 
no use, nor any other, of course. A simple choice is the 11112 model, with four left-handed fermions 
with 0(1) charge qi = 1 and one right-handed fermion with charge = 2. Gauge anomalies cancel, 
because 

Y.iL = Y.iR- ( 51 ) 

i i 

In Fig. 1^, we present a histogram of the change in the phase of the fermion determinant in a lattice 
of L/h = 4 and several values of f/b = 1,1/4,1/8, under a random ^-lattice gauge transformation, 
starting from the trivial configuration. As expected, the change in the phase of the fermion determinant 
gets smaller with the ratio f /b. For arbitrarily small //6, it is clear that the FNN conditions are still 
satished in the full theory. A full unquenched simulation would be needed to determine how small this 
ratio has to be in practice. 

Fermion loops are also expected to induce renormalization of the Yukawa couplings y and r at higher 
orders. For the same reason as before, it can be shown that the corrections to these couplings are 
down by powers of (//6)^. Thus, for small enough //6, no tunning of the bare Yukawa couplings is 
needed to remain in the chiral phase. 

7 Conclusions 

We have presented numerical evidence in favour of the existence of a truly chiral phase in a two-cutoff 
formulation of a quenched 2D Wilson-Yukawa model with a global U{1)l x The Yukawa phase 

diagram is different from the one found in the one-cutoff formulation of this model in that there is 
a new phase for r = 0(1) and ?/ -C 1, where doublers are split from the light spectrum, which is 
composed of massless fermions with chiral quantum numbers. When one of the chiral symmetries is 
gauged, the continuum limit of this lattice model is expected to give rise to a chiral gauge theory. 

The main new feature of the two-cutoff formulation is that the high Fourier modes of the scalar helds 
are naturally cut off at a scale smaller than the momenta of the fermion doublers. As a consequence 

®As explained in the fermion determinant is regulated in a way in which the breaking of gauge invariance is 
restricted to the phase of the determinant. The absolute value can be defined in a gauge-invariant way, otherwise 
one-loop subtractions would be needed. 
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the Wilson-Yukawa term is truly a perturbative coupling for the ligh-fermion mode (p —> 0), while it 
is effectively a large 0{r) Yukawa coupling for the doubler modes {p —>■ tt). Near p —0, the chiral 
fermion propagators are well reproduced by weak Yukawa perturbation theory and have poles at p —> 0, 
while at large momenta they are better described by the strong Yukawa expansion, which predicts no 
doubler poles (since all fermions are massive in this approximation). 

All our results have been obtained in the symmetric or VX phase (PM in 4D) of the scalar dynamics. 
This is the relevant phase for describing a chiral gauge theory. In the case of a spontaneously broken 
gauge theory in 4D, the relevant phase would be a ferromagnetic (FM) one. If a FM counterpart of 
the chiral phase found here exists, it would lead to a successful discretization of the Standard Model. 
This interesting case will be considered elsewhere. 
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Appendix A. 

We present here the explicit formulae for the compact U{1) interpolation in 2D. For details on how to 
construct this interpolation the reader is referred to 0. The methods are similar to those used in 
The Ylattice sites are denoted by s. For an arbitrary gauge conhguration on the 6-lattice, we can 
uniquely write each lattice link variable in the form 

f/^(s) = \A,{s)\<7r, (A. 1) 

where we are neglecting the measure-zero set of lattice helds, where at least one of the link variables 
equals exactly —1. This assignment dehnes an obvious logarithm function for group elements different 
from —1: 

A^ = -i\ogU^, |A^| < TT. (A. 2) 

The link variables on the /-lattice can also be writen as 

= |a^| < TT. (A. 3) 

Dehning the /-lattice sites within each 6-lattice plaquette as x = s -1-^22, where 0 < fi, ^2 < 1, 
the hnal result is, 

Oi = ^ { (1 — ^ 2 ) (Ai(s) — 27rA2(s)) -|- G (Ai(s -1- 2) — 2tiN2{s + 2)) } 

02 = ^ { (1 — ti) (^ 2 ( 5 ) -l- 2nNi{s)) + ti (^ 2(5 -|- i) -|- 2'kNi{s + 1)) } (A. 4) 

6 

Notice that the helds at the boundary of the plaquette satisfy a^ = f /b A^. 

Except for the contribution of the helds, this is simply a linear interpolation, as naively expected. 
The helds are integers dehned by the following equation: 

^ N^{s + (i)-N^{s) = N{s) (A. 5) 
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Figure 16: “Gauge” choice used in solving ([A. 5| ) for compact QED 2 . The corresponding to links 
that are not intersected by the path are zero. 


where N{s) are the winding numbers of the held dehned on the boundary of the plaquette (to be 
more precise the winding of the / —> 0 limit of on the boundary), i.e. 


I[ 


^2('S) + ^ 1(5 + 2) — ^ 2(5 T i) — ^i(s)- 

2^r ■ 


(A. 6) 


with I[y] = nearest integer to y. The existence of these helds is of course related to the compactness of 
the group. In particular, the geometrical dehnition of topological charge is simply given by, 

In this paper we are only concerned with topologically trivial conhgurations since g = 0. Even in 
the case of trivial topology, a non-measure zero set of lattice conhgurations will have N{s) 7 ^ 0. The 
number of non-zero windings grows with the 6 -lattice volume. 

Equation (|A. 5|) has a unique solution in the “gauge” depicted in Fig. All the links that are not 
crossed by the path have the associated iV^(s) set to zero (remember that y is orthogonal to the link 
in our notation). In order to ensure the discrete lattice symmetries, however, we have to also consider 
the interpolations obtained in the maximal gauges that are obtained from that in Fig. 1 by lattice 
symmetry transformations (i.e. 90° rotations and temporal and spatial inversions). There are eight of 
these maximal trees in total. Ideally, the maximal path could be chosen randomly for each 6 -lattice 
conhguration. This is however more complicated to implement, so we have used the hrst method. 

We also need the interpolation of 6 -lattice gauge transformation, dehned as 


u[U^]{x) =u^[U]{x), 


but in the simpler global limit, i.e. 


u[I^]{x) = u‘^[I]{x) = /“'(x). 


(A. 7) 


(A. 8 ) 


where / is the identity conhguration, i.e. G^(s) = 1 or M^(a:) = 1. For the prescription ( |A. 4|) , there is 
a unique solution for to. Dehning 


the solution is. 


G(s) = u;{x) = 

0(^1) G) = ‘^’(s) + (^l('S) — 27iN2{s)) ti + (v42(s) -|- 27iNi{s)) t2 
-|-(A2(s -t- i) -|- 27riVi(s + 1) — ^ 2 ( 5 ) — 27riVi(s)) tit2 


(A. 9) 


(A. 10) 
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with 


(A. 11) 


=-ilog(ll^(s)ll(s + /i)), |A^|<7r, 

and the helds are dehned by (|A. 5|) , (|A. 6|) and (|A. 11| ). 

It is easy to see that 0(0, 0) = <h(s) and a; is thns a smooth interpolation of the gange transformation 
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